Non-factorizable effects in the B — B mixing 
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We analyse the corrections to factorization for the B — B mixing amplitude assuming that 
non-factorizable soft gluon exchanges can be described by the local gluon condensate. Within this 
approximation the (a 3 GG)-correction to the mixing amplitude can be expressed through the B- 
meson transition form factors at zero momentum transfer. The correction is negative independent 
of the particular values of the form factors. 

As a next step, we study these form factors making use of the relativistic dispersion approach based 
on the constituent quark picture. We obtain spectral representations for the form factors in terms of 
the B- meson soft wave function and the propagator of the color-octet bq system. The B- meson wave 
function has been determined previously from weak exclusive B decays, so the main uncertainty in 
our estimates comes from the unknown details of the propagator in the nonperturbative region. We 
obtain AB Bd (m b ) = -0.06 ± 0.035 and AB Bs (m b ) = -0.05 ± 0.03. 

I. INTRODUCTION 

The study of the oscillations in the system of neutral B mesons provides important information on the pattern 
of the CP violation in the Standard model and its extentions (a detailed discussion can be found in Q). The main 
uncertainty in the theoretical description of the process arises from the nonperturbative long-distance contributions 
to the mixing amplitude. In the B — B mixing there are two kinds of such contributions: first, effects related to 
the presence of the i?-mesons in the initial and final states, and, second, corrections to the weak 4-quark amplitude 
due to the soft gluon exchanges between quarks of the initial B and the final B meson. By neglecting the second 
contribution the amplitude factorizes. In this case all nonperturbative effects connected with the meson formation 
are reduced to only one quantity - the leptonic decay constant fg . 

Soft- gluon corrections to the weak 4-quark amplitude give rise to non-factorizable contributions, such that the 
influence of the B-meson structure is no longer described by the decay constant only. The understanding of the actual 
size of the non-factorizable effects and their adequate description in the heavy quark systems is an important and 
challenging task. 

The theoretical analysis of B — B mixing was first performed within the QCD sum rules ^-Q] and later using lattice 
QCD, see ||[| and refs therein. The results from the sum rules are in all cases well compatible with factorization 
(vacuum saturation), within rather large errors. Recent lattice analyses definitely reported negative non-factorizable 
effects, around 5-10% for the central values at the scale n ~ mi,. Still, the errors of the calculations have a similar 
size. Results from these approaches are shown in Table [| 

TABLE I. Comparison of results from various approaches 













Bb 3 1 B Bd 


Ref. 




B Bd (m b ) 


B Bs (mb) 


B Bd 


SR 


20 


1.0 ±0.15 








SR 


ii 


0.95 ±0.1 








Lat 




0.92(4)^^ 


0.91(2)lg 


1.41(6)3 


0.98(3) 


Lat 




0.93(8) 


0.92(6) 


1.38(11) 


0.98(5) 


This work 


0.94 ± 0.035 


0.95 ± 0.03 


1.4 ±0.05 


1.01 ±0.01 



In this letter we analyse corrections to factorization in the B — B mixing due to soft gluon exchanges in the leading 
a s -order, assuming that the main effect of such exchanges can be described by the local gluon condensate [Q. We 
show that within this approximation the parameter ABg which describes correction to factorization, is determined 
by the meson transition form factors at zero momentum transfer. It is important that ABg turns out to be negative, 
independently of the particular values of these form factors. 

As a next step we calculate the relevant form factors. To this end we make use of the relativistic dispersion approach 
based on the constituent quark picture || . We obtain spectral representations for the form factors in terms of the 
B-meson wave function and the propagator of the color-octet 6g-system. 
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Parameters of the model, such as the quark masses and the wave functions, are known from the analysis of the weak 
B-meson form factors [^,0 within the same dispersion approach. For numerical estimates we also need the propagator 
of the color-octet g&-system in the nonperturbative region, so we discuss the possible form of this propagator. The 
lack of reliable information about its details in the nonperturbative region leads to the main uncertainties in our 
numerical estimates. 

The paper is organised as follows: The next Section presents basic formulas for the B — B mixing. In Section III we 
calculate the non-factorizable contribution to the mixing amplitude induced by the local gluon condensate and show 
that the correction is negative. In Section IV we calculate the B-meson transition form factors within the dispersion 
approach. Section V gives numerical estimates. 

II. EFFECTIVE HAMILTONIAN AND THE STRUCTURE OF THE AMPLITUDE 

The effective Hamiltonian which summarizes the perturbative QCD corrections for the B° — B Q mixing has the 
following form (Tl]] 

C 2 M 2 

Ht s B f 2 = ~^f~ (VttVtdf C(n) (bO a d) (bO a d) + h.c. (1) 

where Gf is the Fermi constant, O a = r y a {\ — 75), 7 s = Z7°7 1 7 2 7 3 . The parameter /j, stands for the renormalization 
scale which separates the hard region taken into account perturbatively and the soft region. The Wilson coefficient 
C{p) includes perturbative corrections above the scale \i. The explicit expression for C(fi) can be taken from JlT| . 
The soft contributions are contained in the B-meson matrix elements of the operators in the effective Hamiltonian 

A = (B°\(bO a d)(bO a d)\B°). (2) 

The diagrammatic representation for the amplitude A of Eq. (|^) is shown in Fig |]. 




FIG. 1. Different representations for the B° — B° mixing amplitude. The enclosed quark lines show the order of contraction 
of the spinorial indices. The circles stand for the initial (final) B (B) meson. The second line is obtained from the first line by 
performing the combined color and spinorial Fierz rearrangements. 



Using the language of the hadronic intermediate states one finds that the contribution of the hadronic vacuum to 
this amplitude gives 

(B°\H*r 2 \B°) = l G2FM ^ Ml {V t lV td ) 2 C^)Sl, (3) 

with Jb the leptonic decay constant of the B meson defined according to the relation (0|&7 M 75<i|i?(p)) = i fBPp,- It 
is convenient to parametrize the full amplitude as follows 

(B°\H^\B°) = ~ {v: b V td fC{»)f B B B , (4) 

such that the quantity Bb — 1 = ABb measures contributions of the non-vacuum intermediate hadronic states. 

Using the language of quarks and gluons, the corrections to the amplitude of Fig. 1 due to soft gluon exchanges 
are obtained by inserting the (soft) gluons between the quark lines in Fig 1. Soft gluon exchanges between the quarks 
of the same loop lead to the a s -corrections either to the meson vertices or to the quark propagators. They only 
contribute to the leptonic decay constant fs but do not lead to non-factorizable effects. 

The non-factorizable effects originate from the soft gluon exchanges between quarks of different loops. To lowest 
a s -order these effects are described by the 4 diagrams shown in Fig. 0. 
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b d X -^b a . "^Tb d 

FIG. 2. Soft-gluon exchanges leading to the non-factorizable effects. 



We assume that the main effect of the soft gluon exchange can be described by the local gluon condensate. In the 
next section we demonstrate that in this approximation the (a s GG) -correction to the factorization is negative. 



III. AB IN TERMS OF THE LOCAL GLUON CONDENSATE 

Let us consider the exchange of a soft gluon between different quark loops assuming the dominance of the local 
gluon condensate Jj}. In this case a typical graph of Fig 2 describing the soft-gluon contribution is reduced to the 
product of the three-point diagrams as shown in Fig 0. 



o ^ 




b =( y-mz eu„ ; b 

FIG. 3. Diagrams for T^ {l \ P ,q) (a) and T^ (r) (p,q) (b). 




The corresponding amplitude has the form 



.1'" <r I dqd q 'T?W(p, q )T£?\ P , q ) ^ )4 (2 ^ 4 e-^'-'^(,)i>',(,')), (5) 



where T^a^ (j>, q) with the superscript I stands for the amplitude of Fig [|a with the B in the initial state, and 
T&a' (Pi l') w hh the superscript r for the amplitude of Fig 0b with the B in the final state. 



It is convenient to use the fixed-point gauge x^A h (x) — [p_2|. In this case the gluon potential reads 



A b a (x) = ±x G b af3 + ..., (6) 

where dots stand for terms with higher order derivatives. The average over the hadronic vacuum is performed 
according to the relation 

S bb ' 

( G a0 G a>l3>) = -^r(gaa>gf3f3> ~ g a f3> 9a' ff) (GG) , (7) 

with S aa = N% — 1. Here (GG) is the positive-valued gluon condensate (^GG) = 0.012 GeV 4 §. 
The amplitude of Fig 3 then takes the form 

A W = l J dqdq 'T^'\p^ q )^S( q )^ 

^(G b a0 G b :, pl )^T^( P )^T^, ,( P ), (8) 



where we have denoted 



^T^H P , q )\ q=0 ^^T^l( P ). (9) 

Notice that only the parts of the amplitudes T^' r ) antisymmetric in indices a and (3 give a nonvanishing contribution. 
These antisymmetric parts have the following general Lorentz structure: 
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T aa,f = *P V [fUafjau + if 2 {g a vgp<7 ~ 90v9 a a)] , 

T aa,f = A P V [fUaflav + ifl{g a v9f]a ~ 9/3u9 a a)] , 

rpA(r,q) rpA(l,q) 

aa,/3 act, f) ' 

rpA(r,q) _ rpA(l,q) 

aa,/3 act, [3 ' 

fl = fl fl = -fl (io) 

The superscript q(q) in these expressions denotes the flavour of the quark (antiquark) to which the soft external gluon 
is attached in diagrams of Fig 3. The quantities ff'% are real constants. The last three lines in Eq. ( |l0| ) are the 
consequence of the C-invariance of the strong interaction iS-matrix. As we discuss in the next section, the constants 
/f'f can be represented as specific B- meson transition form factors at zero momentum transfer. 

We have to collect now contributions of the four diagrams in Fig 2. For instance, the contribution to the amplitude 
of the subprocess in which the soft gluons are attached to the b quark in the left loop and the b quark in the right 
loop reads 

T^ b \p)T^$(p)(g aa >gw - g^g*?) = -192M% {{ft) 2 + {fl) 2 ) . (11) 

Similar expressions are easily obtained for other diagrams using the relations ([l^). Taking into account the relevant 
symmetry factors as shown in Fig 1, we finally arrive at the following expressions (we also list the factorizable 
amplitude in the same normalization) 

a, 



A« ~ -8C F M B {^GG) ^ [{ft + fff + {f\ f d 2 ) 2 } , (12) 



with the color factor Cf = ^ N ■ Finally, for ABb wc find the expression 

AB B {^) = [(/* + fff + {f b 2 - fl) 2 ] . (13) 

Jb l j 

Obviously, the correction to the factorizable amplitude due to the local gluon condensate is negative. This agrees with 
all numerical results for AB. Notice that the constants f\ and f 2 depend on the renormalization scale /i. 



IV. FORM FACTORS IN THE DISPERSION APPROACH 

We now proceed to the calcuation of the form factors f\ and f 2 . To this end we make use of the relativistic 
dispersion approach based on the constituent quark picture. Within this approach all observables are given by the 
spectral representations in terms of the soft wave function of the participating mesons. 

We start with the amplitude T„a ah {p, q). It is diaginal in color indices, so we write Tal ab {p,q) = -^=Taa{p,q). 

The triangle diagram for Tal{p, q) is shown in Fig. ||a. Throughout this section we shall omit the superscript /. The 
quark structure of the B-meson is described by the vertex iq~i"f5q2G{s) / \J N c , where s is the invariant mass of the qiq 2 
pair |^]. One should set mi = nib, m 2 = for the calculation of f^ h \ and mi = m^, m 2 — nit, for f^\ 




FIG. 4. Diagrams for T aa {p,q) (a) and T CTa ,^(p) = -^ F T< Ta {p,q)\ q =o (b). 
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We need to calculate T aa ^(p) = -^T aa {p, <?)| g =o @- For the T aa (p,q) shown in Fig |k, its derivative T aat p(p) is 
given by the diagram in Fig ||b. The trace corresponding to the Feynman diagram of FigQb reads: 



Saa;j3 = Sp 



i^{m 2 - k 2 )O a {mi + fci)7a(mi + ki)^p{mi + ki) , (14) 



where O a = 7c (1 — 75)- Since this expression is further multiplied by G a p, only its antisymmetric part in a and f3 is 
necessary. The latter reads 



oA 

°cra;/3 



(ml - fci)Sp (i^{m 2 - fc2)7<r(l - 75)7/3( m i - 
= (m\ - kl) [Aep aau (m 2 kl + m\k 2 ) + ^im\k 2 (g a agi3u - gpagav)] ■ (15) 
The corresponding expression for the Feynman amplitude takes the form 

-1 



;3{P) = TTTTT I dk 1 dk 2 5{p-k 1 - k 2 )—^ 2 , -7 

H (27r) 4 « J \m\ — k\Y\m 2 — k 2 ) 



rpA 



2 ^2) 

x ^^(mafcj + mifca) + 4imih%(g a< rgp u ~ gp a g a u)\ ■ (16) 



We define and / s by the relations 



P u fv(v% 



(2^i Jdhdk 2 S(p- kl k 2 ) (m? - fc?) J K _ ^ = / S ( P 2 ). (17) 
After the fc-integration in Eq. ([!(]) one recovers the structure of the amplitude from Eq. ( |l0| ) 

Taa-lliP) = -V \h tpaav + ifl (gaogfiv ~ g^au)] 5 (18) 

with 

fi = (m 2 - mi)f v + mifs, 

h=m 1 (f s -f v ). (19) 

Let us obtain now spectral representations for the form factors ff and ff. We first notice that the expressions ( |T7j ) 
correspond to the triangle Feynman diagram at zero momentum transfer q = 0. This condition means that both 
q 2 = and p 2 — p' 2 , where p' = p — q, so we can write ff (p 2 ) = ff (q 2 = 0,p 2 ,p' 2 = p 2 ), with i — v,s. 

It is convenient to start with the case q 2 ^ and p 2 ^ p' 2 . Then the form factors fi can be written in the form of 
the double spectral representation [||: 

ff (q 2 ,p 2 ,p> 2 ) = / d \ ^ l2 ^{s,S>, q % (20) 

J s — p z s — p 

where the double spectral density can be calculated for any of the form factors from the Feynman integral. 

As the next step, we set q 2 = 0, but still treat p 2 and p' 2 as independent variables. In this case the double spectral 
representations of the form ( p0| ) simplify to the following single spectral representations 

F 2 , 2N _ 1 [ j_ 1 1 l ra ^s + mf-ml + A 1 / 2 ( s ,mf,ml) 



IAp ' pl) =AM 2 j dS s~p 2 s- p> 2 l ° g \ s + m 2 - m 2 - AV3( S , m \M) . 

l07T z J 8 — p z S — p z S 

The representations for the form factors in the form (|2l]) are valid however only in the region of p 2 and p' 2 (far) below 
the threshold. The reason for that is the following: the representations ( p0| ) and ( f2l| ) are based on the Feynman form 
of the quark propagators. This form is however valid only for highly virtual particles, while in the soft region it is 
strongly distorted by nonperturbative effects. In particular, the pole at k 2 — m? in the propagator of a color object, 
like quark and gluon, is absent. Recently, this was confirmed in a lattice study of the gluon propagator Jl3| . 

The modification of the quark propagator in the soft region leads to the change of the spectral representation for 
the form factors in the region of p 2 and p' 2 near and above the bq threshold. 
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Notice that the quantities l/(s — p 2 ) and l/(s' — p' 2 ) in Eq. ( p0[ ) are the propagators of the initial and final bq 
states with virtualities s and s', and the squared masses p 2 and p , respectively. The nonperturbative effects which 
modify the quark and gluon propagators in the soft region, modify the propagators of the bq states as well. 

We know the character of these changes in the p 2 -channel corresponding to the S-meson: As discussed in H , soft 
interaction of quarks effectively replace the factor l/(s — M B ) with a regular i?-meson soft wave function <Pb(s) in 
the integrand of (f2l|). The B- meson soft wave function is normalized as follows || 

1 f^j^^M.^),, , m ^2^ 



_ dsfeis) 2s °' q, (s - (m b - m q y) = 1. (22) 



The leptonic decay constant is given by the expression 



^/N~ c (m b + m q ) f X 1 / 2 (s,m 2 ,m 2 ) s ~ (m b ~ m ) 2 
Jb = 5-5 / ds<p B (s) . (23) 



8tt 2 

In the color-octet p' 2 -channel the Feynman propagator l/(s — p 12 ) in Eq ( |2l| ) is replaced by the propagator D(s,p' 2 ) 
which involves proper modifications in the soft region. We therefore find the following representation for the form 
factors Mp' 2 ) = Mq 2 = 0,p 2 = M 2 Bl p 12 ) 



/«(P ) = ^2 / ds0 B (s)D{s,p' 2 ) '- 

i- r i2\ 1 f j i / \ / /2m As + m? -m? + A 1/2 (s,m?,m|)\ 

We know that D(s,p' 2 ) ~ l/(s — p' 2 ) at large values of s — p' 2 , and also that D(s,p' 2 ) is finite at s = p' 2 . However, 
we do not know details of D(s,p' 2 ) in the soft region. Motivated by the discussion in we assume here that the 
nonperturbative effects can be described by the following modification of the propagator 



1 

s - p' 2 + Ml 



D ^P' 2 ) = Z (25) 



where Mq is a mass parameter. In order to guarantee the absence of the pole in D(s,p' 2 ) in the heavy quark limit, 
Mq should scale with the heavy quark mass as follows: Mq = 0(Aqcd mq). So it is convenient to write Mq in the 
form 

M 2 = wm d m b , (26) 

where ma — A.qcd is the constituent mass of the light quark, and w is a parameter of order unity]]] 
In the next section we make use of Eqs. (pit) and (p5|) to analyse f v ,s(p' 2 ) for Bd and B s mesons. 



V. NUMERICAL RESULTS. 



The parametrization for the S-meson wave function and values of the quark masses have been determined from the 
analysis of the exclusive processes in P,|lC| . The B- meson soft wave function can be written in the form 



n ^s 2 -(m 2 q - m 2 ) 2 w{k 2 ) 
ViS) = T2 V^(m q -m,) 2 ^ (2?) 

with k 2 — A(s, m 2 , m|)/4s. As found in [1(J a simple Gaussian parametrization w(k 2 )~ exp(— k 2 /2f3 2 ) gives a good 
description of the i?-meson properties. The numerical parameters of the model from flo[| are listed in Table ||. 



1 Assuming that for the color-octet light qq system Eq ( Pq ) remains valid with Mq = wm d , we find D(k 2 = 0) ~ 15 GeV -2 
for w = 1. This is close to the value of the gluon propagator D(k 2 = 0) ~ 18 GeV -2 as found in |l3). This agreement seems to 
be reasonable since one can expect the propagator of the light qq color-octet system to have a structure in the nonperturbative 
region similar to the structure of the gluon propagator. 



G 



TABLE II. Constituent quark masses, slope parameters of the Gaussian wave function, and the corresponding calculated 
leptonic decay constants in GeV from |JL0| . 



m u 


m s 


m b 


Pb 


Pb s 


fa 


f Bs 


0.23 


0.35 


4.85 


0.54 


0.56 


0.18 


0.20 



The procedure of Ref. |9|,|l0[ determines the wave function by fitting the lattice data on the weak transition form 
factors for large momentum transfers at the normalization point fi ~ 5 GeV. Therefore, the i?-meson soft wave 
function and the form factors f s ,v(p' 2 ) at this scale are determined. 

Whereas the i3-meson wave function is known quite well, good information about the details of D(s,p' 2 ) is lacking. 
We therefore use the simple Ansatz ( p5|) for D(s,p' 2 ) in the full range of s and p' 2 and treat w as a free parameter of 
order unity. We expect that the variation of w in the interval w = 0.5 ± 2 provides reasonable error estimates for the 
form factors and AB. Notice that the value of AB corresponding to w = 1 agrees favourably with the recent lattice 
estimates, see Table |[ Fig |j] shows the form factor fg which gives the main contribution to ABb- 




o<^ — 1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 — L 

10 20 30 



FIG. 5. The form factor f s {p' 2 ) for the B d meson vs p' 2 . Solid curves are results of the calculation via the Eqs. ( |24| ) and 
( pB] ) for various values of the parameter w: upper curve - w — 0.5, middle curve - w — 1.0, lower curve - w = 2. The dash-dotted 
curve in the region p' 2 < (mb ± vrid) 2 corresponds to w = 0. The vertical dotted line corresponds to p' 2 — M%. 



The calculated form factors are shown in Table III. The relative magnitudes of these form factors can be easily 
understood taking into account their scaling propertiers in the heavy quark limit: /i^ ~ m b 1 ^ 2 , /j 



-3/2 



1. 



The value of AB is obtained from Eq. @ taking (^-GG) = 0.012 GeV 4 from @. Using the Wilson coefficent at 
the scale fj, = nib we obtain the renorm-invariant factors Bb a and Bb s (for definition see ||) listed in Table |l| 

TABLE III. Form factors and ABB q with q = d for Bd, and q = s for B a . The factors /i,2(p' 2 ) are calculated at p' 2 — M% 
for Bd and at p' 2 = M% for B s . The errors correspond to the interval w=0. 5-^-2. 





B d 


B s 


Ab) _ Ag) 

Jv Jv 


9±3 x 10" 3 


8±3 x 10" 3 




1.1 ±0.2 x 10~ 2 


9±3 x 10" 3 




1.4 ±0.5 x 10" 1 


1.0 ±0.4 x 10" 1 


/r 


7 ±1.5 X 10" 3 


8±3 x 10" 3 


f(b) 
J2 


5.5 ± 1.5 x 10 -3 


5.5 ±2.0 x 10~ 3 




7.5 ±2.5 x 10~ 2 


7.5 ±2.5 x 10~ 2 


Ad) 
J2 


3.5 ± 1.0 x 10" 2 


3.0 ± 1.5 x 10" 2 


AB Bq 


-0.06 ± 0.035 


-0.05 ±0.03 
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Our results are in good agreement with the lattice results, except for a slightly different estimate of the SU(3) 
violating effects in the Bd — Bd and B s — B s cases. This nice agreement gives support to our Ansatz for the description 
of the nonperturbative effects in the propagator of the color-octet qq system. 

We would like to notice that taking the Ansatz (|2^) and ( p6|) for the propagator of the sq system and setting 
Mq = wm s mq, with w ~ 1 allows a parameter- free estimate of the K° — K° mixing. In this case the result is 
much more stable with respect to the particular value of w and therefore to the details of the propagator in the 
nonperturbative region. We obtained this way ABk(1 GeV) = —0.21 ± 0.04 p5| . 



VI. CONCLUSION. 



We considered corrections to factorization in the BB mixing amplitude due to soft-gluon exchanges, assuming that 
the main effect of such exchanges can be described in terms of the local gluon condensate. 

1. It was demonstrated that within this approximation correction of the order a s (GG) to the factorization is negative. 
It can be expressed through the specific B-meson transition form factors at zero momentum transfer. 

2. A relativistic dispersion approach based on the constituent quark picture has been used to calculate these form 
factors. Spectral representations for the form factors in terms of the B-meson soft wave function and the propagator 
of the color-octet qb system were obtained. The behaviour of this propagator in the nonperturbative region was 
discussed. The obtained numerical estimates for ABb favourably compare with the recent lattice calculations. 

Let us point out that the proposed approach can be extended to the description of more complicated problems. In 
particular, it can be applied to the analysis of the non-factorizable effects in non-leptonic B decays. 
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